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Abstract 

We investigate the effects induced by the interactions of the Kaluza-Klein (KK) graviton 
with the standard model (SM) particles on the triple Z°-boson production process at the ILC 
in the framework of the large extra dimension (LED) model. We present the dependence 
of the integrated cross sections on the electron-positron colliding energy y/s, and various 
kinematic distributions of final Z° bosons and their subscqucntial decay products in both 
the SM and the LED model. We also provide the relationship between the integrated cross 
section and the fundamental scale Ms by taking the number of the extra dimensions (d) as 
3, 4, 5, and 6, respectively. The numerical results show that the LED effect can induce a 
observable relative discrepancy for the integrated cross section (Sled), which can reach the 
value of 13.11% (9.27%) when M s = 3.5 (3.8) TeV and the colliding energy ^ = 1 TeV. 
We find the relative discrepancy of LED effect can even reach few dozen percent in the high 
transverse momentum area or the central rapidity region of the final Z°-bosons and muons. 



PACS: ll.10.Kk, 13.66.Fg, 14.70.Hp 
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I. Introduction 



The large extra dimensions (LED) model, proposed by Arkani-Hamed, Dimopoulos and Dvali 
(ADD) [T], is an attractive extension of the standard model (SM) because of its possible testable 
consequences. In the LED model only graviton can propagate in a D = 4 + d dimensional space 
with d being the dimension number of extra space, while the SM particles exist in the usual 
(3 + l)-dimensional space. The picture of a massless graviton propagating in D-dimensions is 
equal to the scene where numerous massive Kaluza-Klein (KK) gravitons propagate in (3 + 1)- 
dimensions. So we can expect that even though the gravitational interactions in the 4 space-time 
dimensions are suppressed by a factor of 1/Mp, they can be compensated by these numerous 
KK states. Therefore, in either the case of real graviton emission or the case of virtual graviton 
exchange, it is shown that the Plank mass Mp cancels out of the cross section after summing 
over the KK states, and we can obtain an interaction strength comparable to the electroweak 
strength [21 [3]. Up to now, many works on the LED phenomenology at colliders have been 
done, including vector boson pair productions and association productions of vector boson with 
graviton [H EJ E] • 

The precision measurements of the trilinear gauge boson couplings (TGCs) are helpful for 
verification of non-Abelian gauge structure, and the investigation of the quartic gauge boson 
couplings (QGCs) can either confirm the symmetry breaking mechanism or present a direct test 
on the new physics beyond the SM. The vector boson pair production processes were extensively 
studied in the SM for probing the TGCs El EJ QUI EE] • The direct study of QGCs requires 
the investigation of the processes involving at least three external gauge bosons. Recently, the 
triple Z°-boson production in the LED model at the LHC was studied in Ref. j!2|. 

The International Linear Collider (ILC) is proposed with the colliding energy yfs = 200 ~ 
500 GeV which would be upgraded to y/s = 1 TeV, and the integrated luminosity is required 
to be 1000 fb^ 1 in the first phase of operation [TT]. The triple gauge boson productions at the 
ILC are important processes in probing the TGCs and QGCs of electroweak theory, which are 
related to electroweak symmetry breaking mechanism. If their observables coincide with the LED 
predictions on the triple gauge boson production processes, it means the coupling of graviton to 
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gauge bosons in the LED model would be the causation. Therefore, the understanding of the 
LED phenomenology in triple gauge boson production processes at the ILC is necessary. 

In this paper we study the LED effects on the process e + e~ — > Z°Z°Z° at the ILC. The 
paper is organized as follows: In section II we present the calculation descriptions for the process 
with a brief review of the related LED theory. The numerical results and discussions are given 
in section III. In the last section a short summary is given. 

II. Analytic calculations 

In the LED model, the extra dimensions on a torus are compactified to a radius R/2ir. The 
relationship between the usual Planck scale Mp to the fundamental scale M$ is 



2(4vr) 



-d/2 



Ml = -f^M^ (2.1) 

where Mp = l/y/G/y ~ 1.22 x 10 19 GeV and Gn is Newton's constant. In our work we use the 
de Donder gauge for the pure KK-graviton part and the Feynman gauge (£ = 1) for the SM part. 
The Feynman rules for the relevant vertices including KK graviton used in our calculations are 
given below [13] . where we assume that all the momenta flow to the vertices, except that the 
fermionic momenta are set to be along the fermion flow directions. 

• G^ K (k 3 ) - i){ki) - ip(k 2 ) vertex : 

" 4 [Vfa + + 7"(fci + h)" ~ V^i + h~ 2m^)] , (2.2) 



G^ K (k 3 ) - ZP(ki) - Z a (k 2 ) vertex : 



IK 



B^ml + {C^p*^ - C^ aT )k lT k w + ^E^ikx, k 2 ) 



(2.3) 



G^ K (k 4 ) - $(ki) - ^(A&) - ZP(k 3 ) vertex : 



ie^ (-frf p + tV P " %YvH (v f - a/7 5 ), (2.4) 

where G^ K , ip and Z^ represent the fields of graviton, lepton, and Z°-boson, respectively, 
e = gsinOw is the positron electric charge, £ is the SU(2) gauge fixing parameter, vj, ay are 



the vector and axial-vector couplings which are the same as those defined in the SM, and k = 
v / 167rG f 7v = V2/Mp where the reduced Planck scale Mp = Mp/^/8tt. The explicit expressions 
for the tensor coefficients are given as 

B i^aP = X -(rf»^ -^ a rfP -t^tT*), 

E^P r7 (k 1 ,k 2 ) = ^ u (KK + K k Z + k i k Z)-[' l f a K k i+v' /p kZkZ + (fi^v)}. (2.5) 
After summation over KK states the spin-2 KK graviton propagator can be expressed as [?] 

(2.6) 



a + 2 



where 



and 



167T H d / 2 ~ l 

D(s) 



k 2 M s d+2 



vr + 2il{k/^fs) 



(2.7) 



/(A/v^ = P / dy j- — 2- (2.8) 

Jo 1 - y 

The integral I(A/y/s) should be understood that a point y = 1 has been removed from the 
integration path, and we set the ultraviolet cutoff A to be the fundamental scale M$ routinely. 
We denote the process of the triple Z°-boson production at the ILC as 

e + (pi) + e~(p 2 ) -> Z°( P3 ) + Z°( PA ) + Z°(p 5 ). (2.9) 

The SM-like diagrams for the above process are depicted in FigQJa). In the LED model, the 
KK graviton can couple to Z°-pair and fermion-pair. We present the four additional pure LED 
Feynman diagrams in FigfTJb). 

In our calculations the developed FeynArts3.4 package |14j is adopted to generate all the low- 
est order Feynman diagrams and convert them to the corresponding amplitudes. Subsequently, 
the amplitude calculations are mainly implemented by applying modified FormCalc5.4 programs 
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III. Numerical results and discussions 



In our numerical calculations, we use the following set of input parameters |17j : 



m w = 80.385 GeV, m z = 91.1876 GeV, a ew (0) = 1/137.036, 



sin 2 9 W = 1 - m 2 w /m 2 z = 0.222897, m e = 0.510998928 MeV. 



We know that the ATLAS and CMS experiments found several SM Higgs-like events at the 
location of Mh ~ 125 GeV |18j . Recently, ATLAS reported the searching for extra dimensions 
by using diphoton events in = 7 TeV pp collisions [19]. The results provided 95% C.L. lower 
limits on the fundamental Planck scale Ms between 2.27 TeV and 3.53 TeV depending on the 
number of extra dimensions d in the range of 7 to 3. The diphoton and dilepton results from 
CMS set limits on M s in the range of 2.5 - 3.8 TeV as d varies from 7 to 2 at 95% C.L. [20]. 
In this work we take M H = 125 GeV, and set M s = 3.5 TeV (or M s = 3.8 TeV) and d = 3 as 
the representative ADD parameters in case otherwise stated. 

In Refs.[T0l lllj there exist the calculations for the SM one-loop electroweak corrections to 
the e + e~ — > Z°Z°Z° process. We make comparison of our LO numerical results with theirs, 
there we adopt the input parameters equal to those in Ref. [10] and use two different packages in 
order to check the correctness of our LO calculations. In Table [T] we present the results of the 
LO integral cross sections in the SM at the = 500 GeV ILC obtained by using CompHEP- 



Figure 1: The Feynman diagrams for the process e + e — > Z°Z°Z° in the LED model, (a) The 
SM-like diagrams, (b) The additional diagrams with KK graviton exchange. 
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M H 
{GeV) 


osm (fb) at LO 


CompHEP 


FeynArts 


Ref.LlOj 


115 


1.0053(3) 


1.0055(1) 


1.0055(2) 


120 


1.0136(3) 


1.0138(1) 


1.0138(2) 


150 


1.0974(3) 


1.0974(1) 


1.0975(2) 


170 


1.2563(4) 


1.2564(1) 


1.2564(2) 



Table 1: The numerical results of the LO integral cross sections for the process e + e — > 
Z°Z°Z° in the SM at the y/H = 500 GeV ILC by using CompHEP-4.5.1 and Fey- 
nArts3.4/FormCalc5.3 packages, and the LO ctsm by using FeynArts3.3/FormCalc5.3 package 
provided in Table 1 of Ref . [10] . 

4.5.1 [16| and FeynArts3.4/FormCalc5.4 packages separately, and the LO ctsm results by using 
FeynArts3.3/FormCalc5.3 package presented in Table 1 of Ref.|10j. We can see that all the 
corresponding cross sections are in good agreement within the calculation errors. 

In the upper plot of Figj2]we present the numerical results of the LO integrated cross sections 
for the process e + e~ — > Z°Z°Z as functions of the colliding energy \fs in both the SM and 
the LED model. The curves reveal that the integrated cross section increases rapidly when 
y/s < 500 GeV and decreases smoothly when y/s > 600 GeV. Obviously, we can find that 
the curve for the LED model decreases more slowly than that for the SM in the region of 
y 7 ! > 600 GeV. We define the relative discrepancy as Sled = aLE v~° ShI to describe the LED 
effect on the integral cross section for the process e + e~ — > Z°Z°Z° , and plot its distribution 
versus \/s in the nether plot of Figj2j We can read out from Figj2] that the LED relative 
discrepancies at the positions of yfi = 500 GeV, 800 TeV, and 1 TeV are 1.15% (0.83%), 
5.69% (4.07%), and 13.11% (9.27%) for M s = 3.5 (3.8) TeV, separately. The LED relative 
discrepancy Sled goes up as the increment of the colliding energy y/s, and the enhancement 
of the cross section oled — &SM becomes larger and larger when the colliding energy y^s goes 
beyond 600 GeV. 

In analyzing the e + e~ — > Z°Z°Z° event, we classify the final three Z° bosons as the leading 
Z°-boson, the next-to-leading Z°-boson and the next-to-next-to-leading Z°-boson according to 
their transverse momenta, labeled as Z\, Z2 and Z%, respectively. The criterion for final Z°-boson 
clarification is based on the conditions of p^ 1 > p^? > p^? . In Figj3] we present the distributions 
of the transverse momenta of Zi^^-bosons in the SM and the LED model, and the corresponding 
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Figure 2: The integrated cross sections for the process e + e~ — > Z°Z°Z° in both the SM and 
the LED model, and the relative discrepancy due to the LED effect (Sled = aLE asM SM ) ' as 
functions of the colliding energy y/s by taking Ms = 3.5 Te^, 3.8 TeV and <i = 3. 

relative discrepancies, d< gf dff jf 2 ° , 5 LE d(Pt' 2 ' 3 ) in * h e conditions of M s = 3.5 Te^, 3.8 TeF, 
y/s = 800 GeV, 1 TeV^, separately. The curves in Figj3] show that the LED effect enhances 
the differential cross section, particularly in the relative high transverse momentum region. 
We can see that in the conditions of y/s = 800 GeV and Ms = 3.5 (3.8) TeV, the relative 
discrepancies Sled(Pt' 2 ' 3 ) can reach the maximal values of 11.08% (8.27%), 13.54% (9.92%) 
and 8.72% (6.17%) separately, while in the conditions of y/a = 1 TeV and M s = 3.5 (3.8) TeV, 
the maximum values of 5 LE d{Pt) are increased to be 29.75% (21.39%), 41.14% (27.96%), and 
19.94% (14.35%), respectively. The curves for the differential cross sections of p^, 1 in Figsj3]^a,b) 
demonstrate the maxima values of ^led{.Pj^ ) are located at the vicinities of pf 1 = 310 GeV 
and p^ 1 = 410 GeV at the y/s = 800 GeV and 1 TeV ILC respectively, while Fig.7 in Ref.[H] 
shows the relative discrepancy increases monotonously with the increment of pf 1 at the LHC. 

The rapidity distributions of final Z-bosons in the SM and the LED model, dasM,LED/dy Zl ' 2Z , 
and the corresponding LED relative discrepancies, 5LED(y Zl ' 2 ' 3 ), with Ms = 3.5 TeV, 3.8 TeV 
and y/s = 800 GeV, 1 TeV are shown in FigsJU^a-f), respectively. The line-shapes of y Zl and 
y Z2 distributions are similar, while the y z ' A distribution seems to be particular. All the curves 
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Figure 3: The distributions of the transverse momenta of Z°-bosons in the SM and 
the LED model, and the corresponding relative discrepancies, defined as 5led(Pt) = 
(^JP - /^p 1 , wi t h M s = 3.5, 3.8 TeV and d = 3. Plots (a), (b) and (c) are for 

the p^ 1 , pp 2 and p^ 3 distributions at the y^s = 800 GeV ILC, while (d), (e) and (f) are for the 
p^ 1 ' 2 ' 3 distributions at the -y/s = 1 TeV ILC, respectively. 
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Figure 4: The rapidity distributions of Z°-bosons in the SM and the LED model, and the 
corresponding relative discrepancies, defined as Sled(v Z ) = (^pP" ~~ d< dy^ ) / ^ " ' Ms — 
3.5, 3.8 TeV and d = 3. FigsJlJa), (b) and (c) are for the y^ 1 , y^ 2 and y^ 3 distributions at the 
-^/i = 800 GeV ILC, while (d), (e) and (f) are for the y^ 1 ' 2 ' 3 distributions at the -^/i = 1 TeV 
ILC, respectively. 
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for specific Z°-boson, such as Z\ (or Z2, Z3), with different values of M5 and \fs look to be 
generally similar. The plots in the FigJ5] show that the curves bulge apparently in the central 
rapidity region, and the relative discrepancies with M$ = 3.5 (3.8) TeV can reach about 8% 
(6%) and 20% (14%) in the central rapidity regions at the yfs = 800 GeV and y/s = 1 TeV ILC, 
separately. The LED effects intensify the cross section evidently in the central rapidity region 
of \y\ < 0.8 where Sled{u Zi ' 2,3 ) are all beyond 5%, and the relative discrepancies are basically 
stable in this region. 

In Figsj5]^a,b) we present the cross sections in the LED model as functions of Ms for different 
d values at the yfs = 800 GeV and yfs = 800 GeV ILC, separately. The SM cross section which 
is independent of these two LED parameters is depicted as a horizontal line. Obviously, Figj5] 
shows that as the increment of either M$ or d, the cross section in the LED model obviously 
decreases and are getting close to the SM prediction. 
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Figure 5: The cross sections in the LED model as functions of M$ with d = 3,4,5,6. The 
additional horizontal line is for the SM cross sections, (a) y/s = 800 GeV. (b) \fs = 1 TeV . 



In the following we consider the inclusive process of 



e + e- -> 



(3-1) 



where the muons are produced by the subsequential Z° decay of Z° — > We take the 

branch ratio as Br(Z° — > /i + ^~) = 3.366% |20j . Since the kinematic distributions of final /x + 
and are the same, we present the distributions of the transverse momentum of fi in both 



the SM and the LED model, 2Ssjji 



- and da £Ji D , and the corresponding relative discrepancies, 
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$led(Pt)> with d = 3, M s = 3.5 TeV, 3.8 TeV and ^fs = 800 GeV, 1 TeF in FigsEJa) and (b), 
respectively. We can read out from the figures that the 5led(Pt) can reach the maximal values 
of 10.83% and 7.78% for M s = 3.5 TeV and M s = 3.8 TeV, respectively when y/s = 800 GeV, 
while in the case of ^fs = 1 TeV the maximum values of 5led{Pt) increase to 29.16% and 



20.37% for M s = 3.5 TeV M s = 3.8 TeV, respectively. 




100 200 300 50 100 150 200 250 300 350 400 450 
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Figure 6: The distributions of the transverse momentum of final muon in the SM and the 
LED model, and the corresponding relative discrepancies (defined as 5led(Pt) — ( ~ ~~ 

^ L )/^ L ) with M s = 3.5 TeV, 3.8 TeV and d = 3. (a) at the ^fs = 800 GeV ILC. (b) at 

they/s = 1 TeV ILC. 

The rapidity distributions of final muon in the SM and the LED model, an d , and 

the corresponding LED relative discrepancies, SLED(y^), with d = 3, Ms = 3.5 TeV, 3.8 TeF 
and \/s = 800 GeV, 1 Te^ are shown in Figsl7^a,b), separately. The curves for the rapidity and 
the relative discrepancy distributions in the central rapidity regions bulge apparently, and the 
relative discrepancies have obvious peaks. There Sled^) for Ms = 3.5 TeV (3.8 TeV) can 
reach 7.5% (5.4%) and 19.1% (13.4%) at the ^ = 800 GeV and = 1 TeV ILC, respectively. 



IV. Summary 

In this paper we study the impact of KK graviton exchange on the scattering process e + e~ — > 
Z°Z°Z° in the LED model at the ILC. This process is very useful in measuring the quartic gauge- 
boson couplings and probing the existence of extra dimensions. We present the dependence of 
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(a) 



(d) 



Figure 7: The rapidity distributions of muon in the SM and the LED model, and the 
corresponding relative discrepancies (defined as SLEDiu^) = ( "dyt^ ~ d< dy^ ) / ^dy^ ) with 
M s = 3.5 TeV, 3.8 TeV and d = 3. (a) at the = 800 GeV ILC. (b) at the y/s = 1 TeV ILC. 

the cross sections in both the SM and the LED model on the colliding energy y/s, and the 
kinematic distributions of final Z° bosons and their subsequential decay products (muons) at 
the ILC. We find the contribution from the KK graviton exchange enhances the cross section 
evidently when y/s goes up beyond 700 GeV. We provide also the distributions of the transverse 
momentum and rapidity of the final produced muon. The results show that the LED effects 
become more observable in highly ranges or the central region of \y\ < 0.8. We also demonstrate 
the relationship between the cross sections in the LED model and the fundamental Planck scale 
Ms with the extra dimensions d are 3, 4, 5, 6, respectively. 
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